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NESTED  BOUNDS  FOR  SOLUTIONS  OF  DIFFERENTIAL  EQUATIONS 


by 

Christian  Fabry 


ABSTRACT 


For  a  particular  class  of  ordinary  differential  equations,  an 
iterative  procedure  is  described  that  gives  a  sequence  of 
nested  pairs  of  lower  and  upper  bounds  for  the  solution. 

Simple  conditions  are  found  under  which  the  bounds  converge 
to  the  solution.  The  method  is  used  to  study  equations  of 
Lanchester1 s  type^  which  have  applications  to  the  deterministic 

i 

description  of  situations  in  naval  warfares 
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INTRODUCTION 


The  dynamics  of  combat  in  general,  and  situations  in  naval 
warfare  in  particular,  can  be  described  deterministically  by 
means  of  differential  equations  usually  called  Lan Chester 
equations,  Lanchester  being  among  the  first  to  investigate 
such  descriptions  of  combat  [Ref,  1]. 

For  those  equations  to  provide  useful  mathematical  models,  it 
is  desirable  that  information  concerning  the  qualitative 
behaviour  of  the  solutions  can  be  obtained  without  too  much 
effort.  Indeed,  even  if  particular  solutions  can  always  be 
found  by  numerical  integration,  it  is  of  more  fundamental 
interest  to  be  able  to  predict,  for  instance,  the  effect  of 
parameters  and  initial  conditions  on  the  final  result  of  the 
combat  or  on  the  value  of  the  variables  at  a  given  time. 

The  qualitative  behaviour  of  the  solutions  is  studied  here 
through  the  determination  of  successive  lower  and  upper  bounds. 
The  first  steps  of  the  iteration,  obtained  by  hand  computation, 
will  provide  bounds  for  the  variables  (which  might  represent, 
for  instance,  force  levels)  that  depend  explicitly  on  the 
parameters.  It  is  then  easy  to  find,  for  instance,  bounds  for 
the  time  required  for  enemy  or  own  forces  to  fall  below  a  given 
value. 

Besides  its  usefulness  in  a  qualitative  analysis,  the  procedure 
developed  here  can  sometimes  be  adapted  for  implementation  on 
digital  computers  in  order  to  determine,  by  rapidly  converging 
approximations,  the  exact  solution.  This  is  justified  when  the 
successive  linear  differential  equations  can  be  replaced  by 
successive  algebraic  equations ;  an  example  of  such  a  situation  is 
given  below. 
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Although  this  work  has  been  inspired  by  the  consideration  of 
equations  of  Lanchester's  type,  the  results  are  by  no  means 
restricted  to  such  equations,  as  appears  clearly  from  the 
hypotheses  made  on  the  equations.  On  the  other  hand,  it  must  be 
mentioned  that  this  note  is  not  concerned  with  the  elaboration 
of  models.  For  these,  the  reader  is  referred,  for  example,  to 
B.W.  Conolly  [Ref.  2]  and  to  K.M.  Mjelde  and  R.R.  Wiederkehr 
[Ref.  3] j  who  analyse  situations  related  to  the  defence  of 
merchant  shipping  against  submarines.  Although  the  work  of 
K.M.  Mjelde  and  R.R.  Wiederkehr  is  also  concerned  with  the 
study  of  the  qualitative  behaviour  of  the  solutions  of  Lanchester 
equations,  its  emphasis  differs  somewhat  from  our  own,  which 
gives  greater  weight  to  the  mathematical  aspects  of  the  problem. 
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1. 


STATEMENT  OF  THE  PROBLEM  AND  NOTATIONS 


Let  us  consider  a  system  of  differential  equations  in  vector 
notation : 


£  =  Ax  +  Y(t^)  , 

Y  £ E  (the  n-dimensional  euclidian  space),  A  is  an  nx n 
matrix,  jr  is  the  derivative  of  JL  with  respect  to  t. 

Systems  of  type  [1]  will  be  investigated  in  this  work,  with  the 
important  restriction  that  the  elements  a^^  of  A  will  be 
assumed  to  be  non— negative ,  except  those  on  the  main  diagonal, 
which  can  be  of  either  sign;  a^  ^  0  (i  ^  j).  If  a  given 
system  does  not  satisfy  those  conditions,  an  appropriate 
transformation  might  force  it  to  do  so.  A  class  of  particularly 
simple  transformations  of  that  kind  is  considered  in  Appendix  A. 
Another  possibility  would  be  to  include  in  I_(t,^)  the  linear 
terms  ai;)  y  .  having  an  unacceptable  sign.  However,  this  will 
in  general  make  our  method  less  efficient. 


As  vectors  in  will 

introduce  the  following 


X 

< 

z 

if 

yi 

< 

z  . 
1 

X 

X 

if 

yi 

z  . 
1 

X 

> 

X 

if 

yi 

> 

X  . 
1 

X 

X 

if 

yi 

X  . 
1 

compared,  it  is  useful  to 
,  for  x>Z>z£En 

(1  s  i  <;  n) , 

(1  i  =£  n)  , 

(1  ^  i  ^  n) , 

(1  ^  i  ^  n) . 


often  be 
notations 


[1] 


4 


2 


RELATIONS  BETWEEN  SOLUTIONS  OF  DIFFERENTIAL  EQUATIONS 


Theorem  1 

If  a^.  s  0  (i  £  j)  and  if  G(t)  is  any  function  satisfying  the 
relation 

Y[t,x(t)]  <  G(t)  for  t  £  0,  [2] 

the  solutions  of  the  differential  equations 

£  =  A^  +  Y(t,^)  ,  ^(0)  =  a  [3] 

z^  =  Ajz  +  G(t)  ,  z_( 0 )  =  b  >  a  [4] 

are  such  that 

y(t)  <  z(t)  for  t  ^  0.  [5] 

PROOF*  Let  us  consider  the  set  of  values  t  ^  0  for  which 
y.(t)  ^  z.(t)  and  let  S  =  Js^<r  S.  • 

If  S  is  not  void,  it  has  an  infimum  §  >  0  and 
a,  ^(t)  <  _z(t)  for  0  ^  t  <  ^  , 

b*  3±  :  y^s)  =  z±(§) 

Then, 

y±(§)  =  (Ax)i  +  Y±  (?,£)  <  (a£)±  +  Gi(§)  =  z±(§) 

and  a  positive  number  6  exists,  such  that 
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< 


for  §  -  6  s  t  <  | 


y±(t)  -  y±(?) 


z±(t)  -  z±(§) 


t  -  § 


Consequently j 


y±(t)  >  z±(t)  for  §  -  6  £  t  <  §  ; 


the  contradiction  with  [6]  shows  that  S  must  be  void. 

Corollary  1 

If  the  condition  b  >  is  replaced  by  the  less  stringent 
condition  b  ^  a,  a  similar  result  is  obtained^  conclusion  [5] 
being  replaced  by 


x(t)  <  z(t) 


[7] 


for  t  >  0 


PROOF.  If  a.  =  b.  for  some  i(l  ^  i  ^  n), 
xi  '  ’ 


Consequently,  for  a  sufficiently  small  positive  number  .-<e .  ,  the 
relation  y^(t)  <  z^(t)  is  valid  for  0  <  t  ^£^.  More  generally, 
a  positive  number  £  exists  such  that 


2;(t)  <  z(t) 


for  0  <  t  ^  e 


Then,  by  Theorem  1, 


^(t)  <  z(t) 


for  t  S  £ 


and  Corollary  1  is  proved. 
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Corollary  2 


If  condition  [2]  is  replaced  by 

Y[t*y(t)]  ^  G(t)  for  t  ^  0  ,  [8] 

conclusion  [7]  must  be  replaced  by 

^r(t)  ^  jz(t)  for  t  ^  0  •  [9] 

PROOF.  Let  £  be  a  monotonically-decreasing  sequence  tending 
n  (  \ 
toward  zero«  Because  of  Corollary  1>  the  solution  z/  '(t)  of 

the  equation 

=  A  z^n^  +  G(t)  +  £  n  h  *  z^n^(0)  =  a.  ^ 

where  all  the  components  of  h  are  equal  to  1^  is  such  that 

_y(t)  ^  z^n^(t)  for  t  ^  0  . 

By  the  continuity  of  the  solution  with  respect  to  a  parameter ^ 
the  relation 

^(t)  £  z(t)  for  t  ^  0 

is  established  for  the  solution  .z(t)  of 
z  =  A z  +  G(t)  ,  z(0)  =  a, 
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If  lower  bounds  are  considered  instead  of  upper  bounds,  similar 
results  are  obtained.  For  instance,  the  following  theorem  can 
be  proved. 

Theorem  2 

F(t)  is  any  function  satisfying 

for  t  s  0  , 

the  solutions  of  the  differential  equations 

±  =  Az  + 

x  =  Ax  +  _F  ( t ) 
are  such  that 

x(t)  £  ^(t)  for  t  2;  0 


Z(°)  “  a  > 

x(0)  =  a  , 


If  a ^  0  (i  £  j)  and  if 
the  relation 

F(t)  ^  Y[t,y(t)] 
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3. 


NESTED  BOUNDS 


Let  us  assume  that^  for  t  ^  Tla  priori11  lower  and  upper 

bounds  and  z^^(t)  are  known  for  the  solution  y(t) 

of  [3];  those  bounds  could  be  obtained^  for  instance,  from 
simple  inspection  of  the  equation.  Assuming  that  Y(t,j^)  is 
continuous,  it  is  possible  to  deduce,  from  and 

z^^(t),  lower  and  upper  bounds  F^^(t)  and  G^^(t)  for  the 


Y[t,^(t)]: 

,(0>(t) 

=  inf 

x(0,(t)  S£(t)  i 

z<°>(t)  lCt’z(t)] 

[10] 

,(0)(t) 

VI 

&  -p 

VI 

[p 

0 

wxi 

II 

,«»(*)  lCt,Z(t)] 

[11] 

of  Corollary  2  and  Theorem 

2,  the  solutions 

J(t) 

of  the  equations 

:(D  = 

A  +  F*-°)  (t)  , 

x(l)(0)  =  a 

[12] 

(D  = 

A  zU)  +  (t)  , 

z^to)  =  a 

[13] 

are  new  bounds  for  y(t)  *  i.e« 
x(l)(t)  *  z(t)  ^  Z(l) (t) 


for  t  ^  0 


Moreover,  if 

.(°)  ^  A  x(o) 

+  F^°kt) 

.  2£(0)(0) 

.(0)  (0) 
z  2:  A  z 

+  G(0)(t) 

,  z(0)(0) 

[14] 

[15] 
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the  new  bounds  will  lie  between  the  "a  priori"  ones: 

x^°^(t)  £  x^^t)  ^  ^(t)  ^  z^^t)  ^z^°)(t) 


This  results  from  the  application  of  the  type  of  argument  used 
in  Theorem  1  and  Corollaries  1  and  2  to  equations  [12]  and  [14] 
and  to  equations  [13]  and  [15]. 


It  must  be  noticed  that,  with  the  conditions  [14]  and  [15],  it 
is  no  longer  necessary  to  assume  that  x(0)(t)  and  z(0)(t) 
are  bounds  for  ^(t),  as  this  results  from  [16],  Because  of 
the  relations  [16],  new  bounds  can  be  found  for  Y[t,^(t)]  that 
lie  between  F^^(t)  and  G^^(t): 


F(1)(t)  = 


(1),  ,  ln5  (11  , 

2L  ^  (t)  s  2^)  2  £  (t) 


G(1)(t)  = 


sup 


(1)  (3)  Y[t,£(t)] 

xUj(t)  *  z(t)  £  z(1)(t) 


(  2  )  (  2  ) 

Then  the  solutions  xv  '(t)  and  zr  '(t)  of  the  equations 


.(2) 
xv  ' 


A  x(2^  +F^1^(t)  ,  X(2^(0)  =  a  , 


z(2)  =  A  z(2)  +  G(1)(t)  ,  z(2)(0)  =  a  , 


are  new  bounds  for  ^(t)  and  satisfy  the  inequalities 
x^(t)  £  x^(t)  ^  y;(t)  ^  z^2^(t)  s  z^^(t) 


because 


x^  ^  A  x^ 2  ^  (t)  , 

£<2>  *  A  z<2>  +G<°>(t)  . 


[16] 
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An  iterative  procedure  is  thus  initiated  which,  at  each  step, 
gives  new  bounds  for  the  solution  y;(t)  of  equation  [3]  j  the 
new  upper  (lower)  bound  being  smaller  (larger)  than  the  previous 
upper  (lower)  bound.  This  procedure  is  explicitly  defined  by 
the  equations : 


*<")  . 


•  (n) 

z  = 


(n) 

+  F^n  1')  (t)  , 

x( n )  (0)  =  _a  , 

[17] 

(n) 

+  G(n  “  1^(t)  , 

z/n^(0)  =  a.  , 

[18] 

together  with  the  equations 


F  (n  1^(t)  =  (  1  i  lnf  /  Y[t,£(t)] 

x(n"  l}(t)  *Z(t)  ^z(n-1}(t) 


[19] 


GU  ~  1}(t) 


,  n  SUP  (n  D  l[t,x(t)] 

^(n-l)(t)  £  z(t)  s  zU"1)(t) 


[20] 


(n  —  1,  •••«) 

If  the  definition  of  F 1  n  ”  ^  (  t )  is  changed  to 


F(n_1)(t)  = 


inf 

x(n"  1}(t)  ^  Z(t)  ^  Z(n)(t) 


Y[t,£(t) ] 


z^n^(t)  being  determined  before  F^n  ~^(t)  ,  a  modified  procedure 
is  obtained  that  will  converge  more  rapidly  than  the  initial  one. 

It  must  also  be  noticed  that,  as  long  as  x^  ^(b)  ^  ^(t)  s  tS  ^(t), 
it  is  possible  to  determine  a  sequence  of  nested  bounds,  even  if 
the  conditions  [14]  and  [15]  are  not  satisfied.  For  that  purpose, 
it  is  sufficient  to  replace  definitions  [19]  and  [20]  by 


11 


IX? 


and 


p(n  -  l)(t) 

^(n  -  1) (t) 


inf 

S2(t) 

sup 

X(n  '  l}(t)  2  j(t)  ^£^n-1^(t) 

(n  =  1,  2 


(n"1^(t)  =  max  [x^n  1^(t),  x(n“2)(t)] 

I(n_l)(t)  =  min  [z^n~1^(t),  z<'n“2^(t)] 

(n  =  2, 

x<°>(t)  =  x<°)(t)  ,  S<°>(t)  = 


•  *  •  •  ) 

9 

•  tit) 
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4* 


CONVERGENT  SUCCESSIVE  BOUNDS 


The  nested  bounds  of  the  previous  paragraph  need  not  necessarily 
converge  to  the  solution  of  equation  [3]*  However,  convergence 
in  this  sense  can  be  ascertained  if  the  function  Y(t^) 
satisfies  a  Lipschitz  condition*  Although  the  notations  in  the 
next  theorem  are  those  of  equations  [17]  to  [20] y  the  result  is 
valid  for  all  the  iterative  procedures  discussed  above* 

Theorem  3 

The  sequences  x^^^)  and  z^n^(t)  5  defined  by  equations 

[17]  to  [20],  converge  to  the  solution  ^(t)  of  [3]  if  a 
constant  L  exists  such  that 


IlK^Zi)  -  K^Zs)  II  s  L  II  Zi  -  z3  II 


f°r  Xu  Z.s  ^  t  S  0  (  ||  x  ||  is  any  convenient  norm  defined 

for  x  £  E  )  . 

—  n 


and  bounded  above, converges  to  a  function  that  will  be  denoted  by 
x^Ct).  Similarly,  the  sequence  z^n^(t)  converges  to  a  function 
z#(t).  Because  Y(t,y;)  is  continuous  in  we  have 


A  x*  +  F*(t) 


x*(0)  =  a  , 


where  F*x-(t) 


X*(t)  sz(t)  s  Z*(t) 


and 


A  z*  +  G*(t) 


z  -5.-  ( 0 )  =  a 


with  G^(t) 


x*-(t)  ^  ^(t)  s  jz*(t) 
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The  difference  f  =  z*  -  x#  satisfies  the  relation 


f  =  A  f  +  G*(t)  -  F*(t)  ,  f(0)  =  0  . 


[21] 


But 5  because  of  the  Lipschitz  condition > 

G*(t)  -  F*(t)  £  L  ||  f  ||  h  , 

where  h  is  a  vector  having  all  components  equal  to  1;  using 
then  Theorem  1  to  compare  equation  [21]  and  equation 


A  =  A  A  +  L  A  h 


A(  0)  =  0  . 


[22] 


it  is  found  that 


f(t)  £  A  (t)  =  0 


However,  by  definition,  f(t)  s  C>  and,  therefore,  f,(t)  =0  oj 


x*(t)  =  z*(t)  =  ^(t) 
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5 


APPLICATION  TO  LANCHESTER-TYPE  EQUATIONS 


5 • 1  Example  Used 

As  an  illustration  of  the  method  discussed  above,  properties  of 
the  solutions  of  the  system 

x  =  -  a  x  -  @  xy  ,  x(0)  =  x0  >  0  ,  [22] 

y  —  -fl'  y -  P 1  x  y  ,  y(0)  =  y0  >  0  ,  [23] 

will  be  studied;  the  parameters  a,  a',  P,  |5'  are  assumed  to 
be  positive. 

This  system  is  a  variant  of  models  proposed  by  B.W.  Conolly 
[Ref.  2]  to  describe  situations  in  naval  warfare.  The  terms 
-Cl  x  and  -a1  y  could  correspond  to  losses  that  either  cannot 
be  attributed  to  the  opposite  side  (such  as  losses  due  to  ship¬ 
wrecks  or  due  to  obsolescence)  or  do  not;  at  least  within  certain 
limits,  depend  on  the  level  of  the  opponent  forces  appearing 
explicitly  in  the  model.  On  the  other  hand,  the  terms  -  3  x  y 
and  -(3*  xy  could  correspond  to  losses  proportional  to  the 
number  of  encounters  between  units  of  the  forces  represented  by 
x  and  units  of  the  forces  represented  by  y. 

The  development  of  the  method  proposed  here  is  in  fact  not 
justified  by  its  application  to  such  a  simple  problem,  but 
rather  by  the  systematic  treatment  it  provides  for  a  class  of 
differential  equations,  of  which  the  system  [22],  [23]  is  only 
a  simple  example.  The  study  of  that  example  will  show  however 
which  type  of  properties  can  be  examined  by  this  method. 
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5. 2  Dependence  on  Parameters  and  Initial  Conditions 


A  conditional  relation  will  be  established  between  x(t)  and 
y(t)  that  gives  some  indication  on  the  effect  of  the  parameters 
and  the  initial  conditions.  In  order  to  establish  that  relation, 
the  transformation  of  variables 


is  used,  the  system  [22],  [23]  being  rewritten  as 


X  = 

-a  X  -  pyG  XY 

x(o)  =  1, 

• 

Y  = 

-  a*  Y  -  P 1  x  0  XY  , 

Y(0)  =  1. 

U  sing 

Theorem 

* 

1,  it  is  found  that 

if 

x(t) 

^  Y(t)  or  £ 

xo 

-  for  t 

a  z  a' 

and  0  y0  s:  0* 

x0 

Similarly,  if  cc  ^  a' '  the  relation 

0  y  ( t)  s  x( t)  C24] 

implies  the  relation 

(By  (t)  >  (3rx(t)  for  t  ^  T 

Possible  values  of  r  for  which  [24]  is  satisfied  could  be 
determined  using  the  bounds  that  will  be  obtained  below* 
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5.3 


Bounds  for  the  Solution 


Bounds  will  now  be  determined  for  the  solution  x(t),  y(t)  of 
the  system  [22],  [23].  It  appears  immediately  from  the 
equation  that 

0  £  x(t)  ^  x0  , 

0  £  y(t)  ^  y 0  ; 

from  this,  it  is  easy  to  find  improved  bounds: 


x0 


-(a  +  P  yoK 

e  ^  x(t) 


-at 

^  X0  e 


[25] 


yo 


-(a'  +  P'x0)t 
e  £  y (t) 


y0 


-a  't 


[26] 


Moreover^  if 


xfl  e 


-(a  +  p  y0)t 


^  X 


1  > 


y0  e 


-(a  1  +  @'x0)t  ^ 


y. 


X, 


/I  "0 

i  .  e .  ,  if  t  ^  miq  (  —  - ,-0 —  In  —  , 
’  va  +  p  y  -  v  ’ 


yox 

0  xl  *  a  »  +  p»  x^  ln 


[27] 


we  have  ^  x(t)  ^  x0  ,  ^  y(t)  ^  y0  and,  consequently, 


x o  e“(a  +  Pyo)t  <:  x(t)  £  x0  e_(a  +  gyi)t  ,  [28] 

-(a1  +P'x0)  /4\  -(a1  +fi'x1)t  Toni 

y0  e  K  0  *  y(t)  s  y0  e  p  l/  .  [29]. 

New  bounds  have  thus  been  found  that  are  valid  as  long  as  t 
satisfies  condition  [27] 
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5.4 


Bounds  for  the  Time  Required  to  Reach  a  Given  Level 


When  studying  Lanchester  equations,  it  is  important  to  estimate 
the  time  required  for  some  variables  to  fall  below  a  given  value. 
Bounds  for  that  time  are  determined  here  in  relation  to  the 
example  considered.  For  instant,  if  T  is  defined  by 
x(t)  =  x  ,  it  results  from  the  inequalities  [28]  that 


1 

ft  +  j9y0 


x, 


in  —  STS 


ft  +  3  yx 


provided  that  yz  satisfies  the  relation 


a’  +  B '  xc 


in 

yi 


a 


i 


Zn 


[30] 


In  other  words,  if  [30]  is  satisfied,  x  reaches  the  level  x1 
before  y  reaches  the  level  y l  .  On  the  contrary,  y  reaches 
before  x  reaches  Xj  if 


i  ,  y0  _ i  ,  x0 

ft  1  +  s’  ln  yx  s  ft  +  e  y0  ln 

Consequently,  the  value  y(xx)  taken  by  y  when  x  equals 
Xj  is  such  that 


ft  1  +  B1  x1 

ft  +  0  y0 


^  'tn 


y(xt ) 


<LLLL£a  . 

ft  +  B y (xj ) 


in 
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5.5 


Numerical  Determination  of  Successive  Bounds 


If  y 0  ^  -g-  and  x0  ^  ~  ,  the  bounds  defined  by  [2 5],  [26] 
satisfy  the  conditions  [14]*  [15] •  Starting  with  those  bounds, 
it  is  then  possible  to  generate,  by  iteration,  a  sequence  of 
nested  bounds  that  can  be  defined  in  the  following  way: 


*(0>(t) 

=  x0  e“(Q  +@yo)t  ? 

y(0)(t) 

I 

CD 

o 

II 

+  g1  x  0  )  t 

- 

I 

X 

1 

_  (n-1)  (n-l)  (n)/  x 

P  *  y  >  x(nj(o)  =x0  , 

[31] 

*  (n) 

y 

I 

£ 

I 

..  (n-l)  (n-l)  (n)/ 

P’  x  y  >r  }( o)  =  y o  , 

132.] 

( n  1 ,  2,  3,  ••••)  • 


The  functions  x^n\t),  y^2n^(t)  are  lower  bounds  for  x(t), 
y(t),  whereas  the  functions  x^n  +  ^\t),  y(2n  +  l)(t)  are 
upper  bounds.  The  differential  equations  [31]*  [32],  can  easily 

be  transformed  into  algebraic  equations.  Indeed,  let  us  assume 
that 


x 


(n~  1)(t) 


N  1  ( n  ,.(n-l). 

E  a.{n-1)  e"Xi  13 

i  =1  1 


i  =  1  1 


or  that 


(n  -  1)  / ,  •w 
x  (t)  y 


(n  -  1 )  /  n  _  v' 


(t)  =  s 

3  =1 


c.U-1) 
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where  the  coefficients  c_.^n  ^  and  exponents  ^ 

easily  be  deduced  from  a_^n  \ /n  ~  \j, ^ 


can 


Then,  if 


Vj  n_1^  ^  a  ,  a'  ( Vj-  )  ,  we  have 


(n) 


N(n)-1  c.(n-1} 


x'  '  (t)  =  £  p  — J 


3=1  V  '"a 


(n  -  l)/t\ 

nr^rr-  •  j  +  *(n)  «"at 


(„1  N(n)-1  o.(n-1) 

*  '<‘>  =  .2  s1  (i-D 

j  =  l  v  .  -  a  ' 


-\j  “  ■*■)  c-h') 

e  J  +  B(n)  e-®'* 


where  N 


B 


(n) 

=  [N(n 

-1’]2  + 

l 

y 

(n) 

N(n)-1 

c.(n-1} 

z 

j  =1 

3 

3 

xo 

(  n  -  1 ) 

v-  -a 

J 

y 

(n) 

IT 

I 

H 

O 

tT 

i 

H 

=  yo  - 

£ 

j  =  i 

0’ 

a 

(n  -  1)  ■ 

Vn-  -a’ 

• 

The  above  described  procedure  has  been  applied  to  the  system 

x  =  -0.2  x  -  0,1  xy ,  x  ( 0 )  =  1  9 

y  =  -0.2  X  -  0.2  xy,  y(o)  =  1 

Starting  with  the  lower  bounds 


*<°>(t)  =e-°-3t,  y<°)(t)  =  e-°‘4t  , 
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three  iterations  have  been  performed,  the  results  of  which  are 
represented  in  Figs.  1  and  2.  The  following  relations  are  shown 
to  hold  between  the  approximate  solutions  obtained  by  the  second 


and  third  iterations: 

lim 

0.69 

tr> 00 

X(3,(t) 

0.72 

1 

lim 

y<2>(t) 

0.38 

t->“ 

0.45 

• 

On  the  other  hand,  it  can  be  proved  that,  for  O^t^lO,  the 
difference  between  x^(t)  and  the  exact  solution  x(t)  is 
less  than  3x1°  »  whereas  the  difference  between  yv  7(t)  and 

_3 

the  exact  solution  y(t)  is  less  than  6  x 10 
The  same  method  has  been  applied  to  the  system 

x  =  0.001  x  -  0.005  xy,  x(0)  =  1  , 

y  =  0.05  y  -  0.02  xy,  y(0)  =  1  , 

where  the  terms  +  0.001  x  and  +0.05y  could,  in  biological 
models,  correspond  to  birth  processes.  The  iterations  have  been 
initiated  with  the  lower  bounds  x  =  0,  y  =  0  in  which  case  it 
is  not  possible  to  obtain  a  sequence  of  nested  bounds  valid  for 
all  t^0.  However,  considering  for  instance  the  variable  x,  it 

is  shown  that,  after  three  iterations,  an  approximation  is  obtained 

_  q 

whose  error  is  less  than  3  X 10  for  O^t^lO;  the  three 
iterations  for  x  are  represented  in  Fig.  3» 


21 


FIG.  1  THREE  SUCCESSIVE  APPROXIMATIONS  FOR  THE  SOLUTION  x(t)  OF  THE  SYSTEM 

x  =  —  0.2  x  —  0. 1  xy  ,  x(0)  =1  , 
y-  —0.2  x  —  0.2  xy  ,  y(0)  -  1  . 


FIG.  2  THREE  SUCCESSIVE  APPROXIMATIONS  FOR  THE  SOLUTION  y(t)  OF  THE  SYSTEM 

x"  «  —  0.2  x  —  0. 1  xy  ,  x(0)  =■  1  , 
y  =  —  0.2  x  —  0.2  xy  ,  y(0)  “  1  . 
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FIG.  3  THREE  SUCCESSIVE  APPROXIMATIONS  FOR  THE  SOLUTION  x(t!  OF  THE  SYSTEM 
x  —  0.001  x  -  0.005  xy  ,  x(0)  -  1  , 

y  -  0.05 y  -  0.02  xy  ,  y(0)  ”  1  . 
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CONCLUSIONS 


It  has  been  shown  above  that  the  solution  of  a  differential 
equation  could  be  approximated  by  successive  lower  and  upper 
bounds •  This  integration  procedure  could  be,  in  principle , 
applied  to  any  differential  equation,  provided  that  "a  priori11, 
bounds  are  known*  However,  as  with  the  somewhat  similar  Picard 
integration  method,  an  appreciable  number  of  iterations  can  be 
performed  only  in  favourable  cases  (essentially  when  the 
successive  approximations  can  be  expressed  as  sums  of 
exponential  functions) • 

If  our  method  has  rather  limited  possibilities  of  use  for 
integration,  it  offers,  on  the  contrary,  a  widely  applicable 
approach  to  the  qualitative  study  of  the  solutions.  Indeed, 
uniformly  valid  bounds  can  provide  information  concerning  the 
value  of  the  variables  at  any  time  and  also  concerning  their 
asymptotic  behaviour.  Those  bounds,  obtained  algebraically,  will 
also  show  the  effect  of  the  parameters  and  initial  conditions  on 
the  solution. 
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APPENDIX  A 


MATRICES  WITH  POSITIVE  NON-DIAGONAL  ELEMENTS 


An  important  restriction  has  been  put  on  the  elements  of  the 
matrix  A  in  equation  [1];  indeed,  the  non-diagonal  elements 
a.  . ( i  7^  j )  have  been  assumed  to  be  positive.  If  this  is  not  so, 
the  above  described  method  cannot  be  applied;  however,  simple 
transformations  might  exist  that  would  transform  a  given  matrix  B 
into  a  matrix  A  with  positive  non-diagonal  elements.  Simple 
transformations  of  a  particular  type  will  be  considered  hereunder, 
namely  transformations  that  when  applied  to  a  vector  y€En  ,  will 
change  the  sign  of  some  components  of  y,  leaving  the  others 
unchanged.  In  other  words,  we  are  concerned  with  transformation 
matrices  T  =  diag  (a,  ,....,  a  )  where  ft  .  =  ±  1  .  We  want  to 
establish  conditions  under  which  such  a  matrix  T  exists  that 
transforms  a  given  matrix  B  into  a  matrix  A  =  T  B  T  with 
positive  non-diagonal  elements,  i.e.  with  elements  a^.  ^0  (i^=j)  . 

Theorem  A.1  For  a  given  matrix  B,  a  transformation  matrix 
T  =  drag  (a  x ,  ......  ft  ),  (ft  ^  =  ±  1),  exists  that  transforms 

B  into  a  matrix  A  =  T  B  T  with  strictly  positive  non-diagonal 
elements  if,  and  only  if,  the  elements  b^.  of  B  are  such  that 

>  0  for  k  ^  1,  t^l,  k  ^  i  ; 


b.lk 

bk  l 

hll 

b 

b 

>  0 

Dlk 

kl 
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Necessary  Condition  s  By  hypothesis ^ 


anl  a,  o  a. 


*lk  °k l  all 
consequently 


^  0  for  k  ^  1,  £  ^  1,  k  ^  i  ; 


2  2  2 


®1  ak  al  blk  hkl  hll  >  0 


and  t>lk  b^  b^  >  0  for  k  7^  1,  1^1,  k  ^  l 

Similarly, 


aik  aki  >  0 


or 


2  2 


ai  ak  blk  bkl  >  ° 


Sufficient  Condition  t  Let 


«1  =  1  j  ak  =  Sgn  (blk)  (k  ^  1) 


Then, 


lk 


'll 


kl 


sgn 

bik 

> 

0 

sgn 

<*>!*> 

\i 

sgn 

<bu) 

b^i 

> 

0 

sgn 

<bik) 

sgn 

(b 

I*,) 

bk  l 

sgn 

<bik> 

sgn 

(b 

11^ 

bk  l 

(k  ^  1)  ; 


(-f-  t4  i)  ; 
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If  we  want  to  find  a  transformation 
elements  of  A  =  T  B  T  are  positive 
conditions 


T  such  that  the  non-diagonal 
(a±_j  >0,  i  £  3 )  ,  the 


bjk  \l  hlj  *  °  £or  k  t  3’  1  t  k  f  * 


and 


s  0 


are  necessary  and  sufficient. 
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